Field-angle-resolved anisotropy in superconducting CeCoIn 5 using realistic Fermi surfaces 
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We compute the field-angle -resolved specific heat and thermal conductivity using realistic model band struc- 
ture for the heavy-fermion superconductor CeCoIns to identify the gap structure and location of nodes. We use a 
two-band tight-binding parametrization of the band dispersion as input for the self-consistent calculations in the 
quasiclassical formulation of the superconductivity Systematic analysis shows that modest in-plane anisotropy 
in the density of states and Fermi velocity in tetragonal crystals significantly affects the fourfold oscillations 
in thermal quantities, when the magnetic field is rotated in the basal plane. The Fermi surface anisotropy sub- 
stantially shifts the location of the lines in the H-T plane, where the oscillations change sign compared to 
quasicylindrical model calculations. In particular, at high fields, the anisotropy and sign reversal are found even 
for isotropic gaps. Our findings imply that a simultaneous analysis of the specific heat and thermal conductivity, 
with an emphasis on the low energy sector, is needed to restrict potential pairing scenarios in multiband super- 
conductors. We discuss the impact of our results on recent measurements of the Ce-1 15 family, namely CeTIn 5 
with T=Co,Rh,Ir. 

PACS numbers: 74.25.Uv,74.20.Rp,74.25.Bt,74.25.fc 



I. INTRODUCTION 

Many heavy-fermion and other novel superconductors are 
thought to possess nodes in the gap function on the Fermi sur- 
face. Since the gap shape is directly related to the symme- 
try of the pairing interaction, knowing the position of nodes 
can shed light on possible pairing mechanisms. Magnetic 
field-angle-resolved specific heat and thermal conductivity ex- 
periments are able to provide detailed information about the 
anisotropy of quasiparticle excitations near the Fermi surface, 
and hence help identify the nodal directions in the bulkPEl 
To implement this procedure it is necessary to have high- 
precision probes that detect small variations under changes 
of the direction of the applied field. A series of remarkable 
experiments proved already the viability of this approach!^ 
However it has proved non-trivial to interpret these exper- 
iments in general. The oscillations in physical quantities, 
as a function of the field direction, change sign depend- 
ing on the magnitude of the applied magnetic field and the 
temperature The location of these inversion lines de- 
pends sensitively on the topology of the Fermi surface and 
the material- specific details of the (multi-) band structure C^iiiJ 
Obviously, this calls for the development of theoretical tools 
that take material-specific properties into account. Further- 
more, it suggests that a quantitative and unambiguous identi- 
fication of the structure of the superconducting (SC) gap re- 
quires the incorporation of realistic Fermi surface (FS) prop- 
erties. 

The unconventional heavy-fermion superconductor 
CeCoIn 5 is an ideal candidate for testing field-angle-resolved 
probes due to the existence of large high-quality crystals and 
accessible temperature and field ranges. Early field-angle- 
resolved thermal conductivity and specific heat measurements 
were controversial on whether CeCoIns has a superconduct- 
ing gap with d x 2_ y 2 or d xy symmetry.^1121 Recent specific 
heat measurements observed the predicted inversion of the 



oscillations at low temperature. 6 This seemed to have settled 
the dispute in favor of d x i_ y i pairing symmetry. 

In this paper we incorporate first-principles electronic 
structure calculations to obtain the realistic tight-binding 
parametrization for Ce-1 15 (CeTIn 5 with T=Co,Rh,Ir) ma- 
terials that reproduce the Fermi surface (FS) topology and 
yield the Fermi velocities, and the density of states (DOS) 
at the Fermi level. We use this FS parametrization as in- 
put for self-consistent calculations of thermal properties in 
the extended Brandt-Pesch-Tewordt (BPT) approximation of 
the quasiclassical Eilenberger equation .™ Use of the tight- 
binding parametrization allows for a numerically efficient 
computation, while keeping the essential character of the low- 
energy band structure that reflects on the hybridization be- 
tween Ce 4/ and In 5p states. Within this framework, we 
consider candidate s- and d-wave order parameters, perform 
a systematic study of the angle-resolved specific heat coeffi- 
cient, 7 = C/T, and thermal conductivity, k, in a magnetic 
field rotating in the Ce-In basal plane. Finally, we construct a 
field-temperature phase diagram of the fourfold oscillations. 

The main results of our calculations, which are applicable 
to a wide range of systems with tetragonal point group sym- 
metry, are: (1) For isotropic gap (s-wave) we find that mod- 
erate FS anisotropies are sufficient to introduce field-angle- 
dependent oscillations in the specific heat and thermal con- 
ductivity in the superconducting state over a significant range 
of temperatures and at intermediate to high magnetic fields. 
In addition, the inversion of the oscillation pattern as a func- 
tion of temperature shows that oscillations are not simply a di- 
rect consequence of the anisotropy of the upper critical field. 
Therefore not all such oscillations at intermediate fields can 
be taken as proof of strong anisotropy in the superconducting 
gap. This result agrees with our recent numerical study of the 
iron-based superconductor TFe2Se2.^(2) The complex field- 
angle dependence of the specific heat and thermal conductiv- 
ity for systems with anisotropic Fermi surfaces suggests that 
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comparison of both quantities with material- specific theories 
is required to identify the pairing symmetry and gap struc- 
ture. This is already important for materials, where the Fermi 
surface anisotropy is moderate, as is the case for the Ce-115 
family. 

The rest of the paper is arranged as follows. In Sec. II, 
we present our tight-binding representation of the two FSs for 
CeCoIn 5 . Detailed analytical and computational formalism of 
the field-angle resolved specific heat and thermal conductivity 
calculations is given in Sec. III. The results of the temperature 
and magnetic field dependence of these quantities, and their 
relative sign reversal in the field-angle oscillation for s-and 
two d-wave pairing symmetries are given in Sec. IV. Some 
comparison with the available data for CeCoIn 5 , CeRhIn 5 , 
and Celrln 5 is also included. Finally, we conclude in Sec. V. 



II. ELECTRONIC STRUCTURE 

First-principles calculations of CeCoIn 5 demonstrate that 
the bands crossing the Fermi level are dominated by strongly 
hybridized 4/ electrons of the Ce atom with weak overlap 
coming from the bp orbitals of the In atom.LL^^U j n m j s WO rk 
our basic aim is to parameterize the true shape of the FSs only, 
while the overall dispersion feature at higher energy is irrele- 
vant for thermodynamic and transport properties. Therefore, 
we use an effective tight-binding model of the lowest energy 
of three 4/— orbitals in a tetragonal lattice. As we are only 
interested in the eigenvalues and not the eigenvectors of each 
band, we absorb the orbital symmetry of contributing orbitals 
into the tight-binding hopping parameters, which makes all 
bands decoupled from each other. With this motivation we 
write the tight-binding dispersion including up to third nearest 
neighbor hopping on the x — y plane and only nearest neighbor 
hopping along c— axis to obtain 

£ fe = -2 y^jtjCj + t 2i c 2i ) - 4:t xy c x c y - E F . (1) 

i 

Here c ai = cos (aki) with i = x,y, z. E F is the Fermi en- 
ergy. We obtain the values of the tight-binding parameters 
after fitting to first-principles dispersions by Ref. 19 shown 
in Fig. 1(a): (t x =t y ,t z ,t 2x =t2y,t xy ,E F )=(-0A2 r 0.05,0,0.09,- 



0. 55), and (-0.17,0.06,0,0.15,-0.47) in eV for the a and /3 
band, respectively. Note that t 2z = 0. The other two bands 
crossing the Fermi level have small areas and are not further 
considered in our two-band model description of CeCoIn 5 . 

The a and f3 bands give two concentric electron pockets at 
the zone corner (M-point), see Fig. 1(b). The k z dispersion 
of each band is more interesting and needs special attention. 
Along the k z direction, both a and (j FSs are more like corru- 
gated cylinders: the a FS has a narrow waist at k z = 0, while 
the /? FS has a belly. Note that only nearest neighbor hop- 
ping along the c axis is sufficient to obtain the qualitative k z 
dispersion of all bands in agreement with the ab-initio band 
structure^ and dHvA experiments 20 [see Fig. 1(c)]. The op- 
posite sign of the t z parameter is responsible for the opposite 
shape of the a and f3 FSs (narrow waist vs. belly). 

III. THEORY AND COMPUTATIONAL METHOD 

For magnetic field H applied at angle a with respect to 
the (100) direction, we compute the field-angle induced su- 
perconducting DOS per spin, N n (uj\ H) (ban d index n = 

1 , 2) by solving the Eilenberger equation 3 8 9 23 within the ex- 
tended BPT quasiclassical approximation pSHD The BPT ap- 
proximation implies a uniform field H over the unit cell of 
the Abrikosov vortex lattice (unit-cell averaged Green's func- 
tion). This produces quantitatively correct results near the up- 
per critical field, and continues to yield semi-quantitatively 
correct descripti on over the range 0.bH c2 (T) < H < H c2 (T) 
for isotropic gap j 24 * 27 * 2 ^ and to much lower fields for no dal and 
strongly anisotropic gaps in single-band models I322EQI 

Here we summarize the key steps of the calculation, and 
highlight the main technical differences between the single- 
and multi-band systems following Refs. |U |9] and [23] The 
main object of interest, the quasiclassical Green's function, is 
assumed to be diagonal in the band space (n = 1, 2), since 
bands are well separated in the Brillouin Zone, and have the 
4x4 Gor'kov-Nambu matrix structure corresponding to sin- 
glet pairing in each band, 

6=f? ) , g n =( > f %(J2f A. (2) 

The Green's function in each band satisfies the Eilenberger 
equation for given Matsubara frequency iuo v = mT(2v + 1), 
which has a simple commutator form!^ 



+ iv n (k f ) ■ V R g n (R,k f ;iuj u ) = 0, (3) 

where the Fermi velocity in band n is denoted by v n (kf), with the wavevector kf on the respective FS. Since this is a homoge- 
neous equation, it has to be complemented by the normalization condition of the Green's functions: 

g n (R,kf]iuj v ) 2 = -7T 2 . (4) 

Furthermore, the off-diagonal Green's functions and self-energies are related by symmetry f (R,kf;iu u ) = 

f n (R, -k f ; iu u y = f n (R,k f] -iuj v f- A^ p (R,k f ; iuj v ) = A™p(H, -k f ; iu v )* = A™*(R,k f] 

The equations for the Green's functions in two bands are coupled indirectly through the self-energies entering the Eilenberger 
equation. The scattering of quasiparticles off impurities with concentration ni mp is taken into account via the self-energy in each 



(iuj u + -v n (k f ) ■ A(R)) f 3 - A n (fl, k f ) - <5-; mp (^), g n (R,k f ; iuj u ) 
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FIG. 1 . (color online) Electronic structure of CeCoIns : (a) Tight-binding fits of the two most relevant bands a (red) and (3 (green) to the 
electronic dispersions of CeCoIns calculated in the local density approximation by Maehira et al!^(b) a and /3 FSs at three representative k z 
values, colored by corresponding Fermi velocity from low (green) to high (yellow), (c) Three dimensional rendering of the computed FSs for 
a and /3 bands (bottom panel) compared with the dHvA experiments 20 (top panel). The a FS has a narrow waist, while the ft FS has a belly. 
The color map of the calculated FSs gives the anisotropy of the magnitude of the Fermi velocities ranging from low (blue) to high (red). 



band, <r^ mp , which is evaluated in the T-matrix approximation for the two-band system 



D + ?,™ p ia 2 A i Z np 



r H j impt"nn i 



f=U + U(N f (k f )G(kf)) F gt. 



T = 



U = 



The f matrix and the impurity scattering potential have the following structure in band space: 

Un U12 

The angular brackets denote the integral over one or the other Fermi surface, as appropriate, e.g.: 

d 2 k f Nf tn (kf) g n (kf) 



til ti2 
hi h2 



{N f (kf)G(k f ) ) FS = diag„ =1>2 



UFS n 



(5) 



(6) 



(7) 



and the corresponding normal-state DOS at the Fermi level is Nf iTl (kf) ~ l/\v n (kf)\. Sometimes we will omit the subscript 
FS for brevity. 

For each T and H the order parameters are calculated self-consistently from the coupled gap equations of the two-band model 
A n (R,k f )=Tj2J2 ( v nm(k f ,k/) N ftm (k/) f m (R,k/i UJ V )) (8) 



FS 



We use a factorized pairing potential at the Fermi surface as V nm (kf,k/) = V nm y n {4>) y m {<l>')> wim the basis 
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function that depends only on the azimuthal angle, see Figs. [2] 
and [3] This means that the order parameters are also factor- 
ized, Ai^(kf) = Ai^^i, 2 (</>)• We couple the bands, for 
simplicity, by purely interband pairing V\ 2 = V 2 \ = — V, 
and assume same symmetries and angular variations on both 
bands 3^.(0) = 3^2 (</>)• This ensures that the order parameters 
in the two bands are strongly coupled, and the temperature and 
field dependence of both gaps is similar, while keeping the 
number of parameters minimal. While there are indications 
that in CeCoIns there is a small excitation gap that closes at 
very low fields, of order 0.1% of the upper critical fieldp^it 
seems likely that this gap is proximity induced on the parts of 
the Fermi surface with low /-electron content that we do not 
consider here. Since all the experiments measuring the field- 
angle anisotropy are carried out at fields, which are at suffi- 
ciently high H/H C 2, we consider only Fermi surface sheets 
with strong pairing and large gaps. It is worth mentioning that 
the interband pairing captures both nodeless s ± and nodal d 
wave pairing scenarios. 

Generally, for arbitrary interaction matrix V nm , the cou- 
pled gap equations support two solutions for the amplitudes 
(Ai,A 2 ). The physical solution corresponds to the high- 
est transition temperature T c o, that is, the greatest eigenvalue 



Vm ax of the interaction matrix, 



VniNf^y*} v 12 (N ft2 yi) 
v 2 i(N fil y*) v 22 (N f , 2 y 2 ) 



e 2 J "™ y e 2 

(9) 

The effective interaction strength V max and the cutoff Q c can 
be eliminated using standard techniques in favor of the bare 
transition temperature, T c o = 1.13f2 c exp ( — 1 / Vmax ) and 
the gap amplitudes in different bands are given by the eigen- 
vector of the interaction matrix, 



Ax 
A 2 



A. 



(10) 



Upon projecting out this vector from Eq. ([5]), the system of the 
self-consistency equations is reduced to a single equation for 
the order parameter A of the dominant instability. 

Since the only coupling between bands is via the self- 
consistency equations of the order parameter and the self- 
energies, the solutions for the propagators in each band 
can be formally obtained from the transport equation ^ 
with given A n and a n in the same way as for single- 
band systems We express the gradient term via the rais- 
ing and lowering operators {a\a) for the vortex solu- 
tions corresponding to the superposition of different har- 
monic oscillator functions!^ v n (fc/) • (V^- i 2 fA(R)) — 
y=X { v n-(kf)a^ -v ni+ (kf)a) . The (x,y) projections of 
the Fermi velocity on the plane perpendicular to the direction 
of the field z have to be rescaled by the anisotropy factor Sf, 

(kf)y/S~f. (11) 

The relevant parameter that determines the excitations in the 
SC state at a particular point on the Fermi surface is the com- 
ponent of the (rescaled) Fermi velocity normal to the applied 
field, 



On (*/) = Jvn,xQ*f) 2 /Sf + Vy(*/) 2 Sf- (12) 



The corresponding energy scale is 



Vf, n (4>, H) 



<(kf) 



2A 



(13) 



where A = (hc/2\e\H) 1 / 2 is the magnetic length, which is of 
order of the intervortex distance, and </> is the FS angle with 
respect to the k x axis. The anisotropy parameter Sf is chosen 
to give the correct form of the vortex lattice in the linearized 
Ginzburg-Landau (GL) equations for A. This allows us to 
consider only the lowest Landau level, 8 



E, , e ik vV s fy 
ci n) e . = 



'x-\ 2 y/S~fk y 



(14) 

For tetragonal symmetry this parameter depends on the an- 
gle that the applied field makes with the symmetry axis c (in 
this paper Oh = tt/2), 



Sf = 



I cos 2 Oh 



1^ 
K ± 



sin vh • 



(15) 



Here K\ \ (along c-axis) and K± (in-plane) are the coefficients 
of the gradient terms in the GL expansion for the gradients 
along the c-axis and in the afr-plane respectively. For our 
two-band system they depend on the degree of mixing of the 
bands in a particular superconducting state A. For the state 
(Ai,A 2 ) = (eiA,e 2 A) they are determined by the right, 
e = (ei, e 2 ) T , and left, e' = (e' l5 e 2 ), eigenvectors of the in- 
teraction matrix in Eq. ([9]), corresponding to eigenvalue V max 
with e' • e = 1. 



(ei , e' 2 ) 




Via 


Vmax 


\v 21 




(ei , e> 2 ) 


(Vu 




Vmax 


\V 2 i 


V22 



ei (vlN fl y*) 
e2 (v 2 2c N f2 yi) 

ei {v 2 la N fl yl) 
e 2 (v 2 2a N f2 yi) 



,(16) 



.(17) 



With these remarks in mind, we can directly use the single- 
band results for the unit-cell averaged Green's functions in the 
single Landau level approximation [we follow the notation of 
Eqs. (46)-(48) in Ref.|5Q: 



g n (iUv,kf,H) 



j A n A 



= , (18) 



f n (iuJ Vl k f 'H) = ign^W (^f) An. (19) 

Here iuj Ujn = iuo v - Y, % ™ p (iu u ,kf\ H) and A n = A n (fe/) + 
A^ mp (zcj iy , kf, H) are the Matsubara frequency and the order 
parameter renormalized by the impurity self-energies in each 
band n, & % ™ v . W'(z) is the first derivative of the complex- 
valued function W(z) = exp (— z 2 )erfc(— iz). One can fur- 
ther cast this in a form similar to that of a uniform supercon- 
ductor by introducing the new self-energy H n according to 
iy/n/v 2 ^ n W^(iuj u /vf yn ) = (iujj, - H n )" 2 . The effective 
self-energy S n now contains effects from both the impurity 
scattering and the effects of orbital magnetic field. 
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In contrast to the Doppler shift approximation, both the real 
and the imaginary parts of S n contribute to the SC DOS, and 
their interplay as a function of energy, H and T, determine 
the sign reversal in the fourfold oscillation of the SC DOS. 
These effects have been extensively studied earlier using a 
single quasi-cylindrical FS and nodal gap, and a minimal 2D 
model for two-band systems, see for example Refs. 3,8,9] and 



The transport and thermodynamic coefficients are calcu- 
lated by using the retarded Green's functions through ana- 
lytic continuation, iuo v —> uj + iO, in the propagators found 
above. We begin with the total electronic specific heat from 
both bands, C = C\ + C2, which is given by the deriva- 
tive of the net entropy, C = T (dS/dT). Because the low- 



temperature approximation given by^ 

uj 2 (N n (uj,k f ;H)} FS 



C n (a) 



duj- 



4T 2 cosh(uj/2T) 2 



(20) 



remains valid almost up to the normal-state transition region, 
it can be employed to describe the behavior of the heat ca- 
pacity over most of the phase diagram. Detailed numerical 
calculations show that the high-temperature sign reversal line 
is robust, but will be shifted to slightly higher temperatures by 
approximately 0.1T c0 for the FS parametrization considered 
here compared to the calculations using the low-temperature 



approximation in Eq. (20). 

Next, we consider the total electronic thermal conductivity, 
which is the sum of the contributions from both bands, k = 

ftl + K>2, Wit 



duj- 



UJ 



2T 2 cosh(cj/2T) 2 
x (v*(k f ) 2 N n (uj,k f ; H)r n (uj,k f ; H)) FS • 



(21) 



Here the field-induced SC DOS per spin in each band, N n (uj, fc/)/AT/ ?n (fc/) = — im g~\uj, k f ; ri 
for the spin degeneracy, and the transport lifetime is due to both impurity and vortex scattering} 9 ! 23 ! 29 ! 



Im g R (uJ, kf \ H ) n /ir, the factor 2 accounts 



1 



2r n (u),kf]H) 



= -lmY?™v{uj,k f -H) + ^- 



1 Im \g*(u>, k f ; H) W(uj/\v fjn (k f] H) 



Vf in (kf]H) 



\A n (k f ;H)\ 2 (22) 



When T — > we recover the standard expressions for the Sommerfeld coefficient, j n = C n /T — >> |7r 2 (iV n (0, fc/; H))fs> 
and for the linear coefficient of the thermal conductivity /T — >> |7r 2 (^(fc/) 2 Af n (0, fc/; H)r n (0, kf, H)) F s- Since the 
Green's function, given by Eq. (fT8l), takes the standard BCS form at H = 0, we also recover the universal thermal conductivity 
for gaps with nodes on the FS P 5tt39 1 At low fields the approximation breaks down, but for nodal superconductors it provides 



a good interpolation from low to high fields, and, in the regime 1 <C l /r ? ; mT? A n <C H/H C 2 reproduces the well-known y/H 



field-dependence of the density of states in d-wave superconductor^^ up to logarithmic corrections ! 8 I 9 I 29 I 

Since the function x 2 /cosh(x/2) 2 peaks at x ~ 2.5T, the anisotropy of the heat capacity at low temperatures is qualitatively 
determined by the anisotropy in the total SC DOS, N(uj = 2.5T, fc/; H). Using the expansion of the error function, we obtain 
two limiting values for W'(z): W'(0) = 2i/^ and W'(z » 1) w -i/^z 2 . Thus the SC DOS for each band n becomes 



N n (uj;H) = (N n (uj,k f ;H)) FS 




-1/2N 



-1/2 \ 



FS 



UJ > Vf iTl . 



FS 



(23) 



The first line in Eq. ( [23) only makes physical sense when the 
BPT approximation is valid at low energies, i.e., for nodal and 
strongly anisotropic gaps. In that case at low T (low energy) 
and low fields, where A n (fc/;i2") only weakly depends on 
the direction of the field, the SC DOS depends predominantly 
on the orientation of v/ ?n (fc/; H) relative to the minima of 
A n (fc/; H). At uj = the inversion of the SC DOS as a 
function of the field for nodal gaps can be obtained in analogy 
with Refs. [8]and|42l 

At higher energies, the second line of Eq. ( [23] ) has the 
BCS form apart from the replacement of the bare energies 



and gaps by their impurity renormalized counterparts. There- 
fore the field-angle variation enters via the anisotropy of these 
self-energies as well as via the field dependence of the gaps, 
A n (fc/ ; H) , for determining the anisotropy of the upper criti- 
cal field. The latter effect is only relevant in the vicinity of the 
transition to the normal state, where the result is valid for both 
nodal and nodeless gaps, including the fully isotropic situa- 
tion. Crucially, for anisotropic Fermi surfaces the anisotropy 
in the self-energies and the order parameter is weighted by 
the normal-state angle-dependent DOS, 7V/ n (fc/), leading to 
a complex behavior including the switching of the minima and 
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FIG. 2. (color online) Fermi surfaces (FSs) and Fermi velocities at 
k z = (left panel) and k z = tv/c (right panel). The relative magni- 
tude of the Fermi velocities (in arbitrary units) is given by the length 
of the (red and green) arrows in the top panels. The bottom panels 
show the Fermi velocities along the Fermi lines of the corresponding 
k z slice. 



maxima found in this work. However, in this regime the en- 
ergy width of the Fermi weighting factor in the integral ex- 
ceeds the gap amplitude and a full numerical evaluation is 
required. Our results are consistent with the general obser- 
vations based on such an expansion. 

For each pairing symmetry, the coupled order parameters 
are computed self-consistently at each temperature T and for 
a given value of the magnetic field H applied at the angle a to 
the (100) direction. We calculate the field-angle oscillations in 
the H-T phase diagram for a mesh of 35 field points between 
zero and H c2 , 100 temperature points from zero to T c0 , and 31 
field- angle points a from zero to 90° to extract the anisotropic 
terms in the heat capacity and thermal conductivity. For all 
the calculations, we consider purely intraband impurity scat- 
tering, U12 = U21 = 0, t^ii = U22 = uo, in the clean limit, 

o /_ A7 r\ rvi . . c\_rn „.l rn 

cO 



2n imp /7rN fl = 0.01 x 2ttT c0 , where T c 



2I\ — l/ri mPj i 
is the bare transition temperature, Nfi is the density of states 
on the first band (a-band), and the scattering phase shift is 
chosen to be S = arctan(7riV/iMo) = tt/2 (unitarity limit). 



IV. RESULTS 

A. Field-induced superconducting DOS anisotropy and the 
role of Fermi surface topology 

As discussed above, an important aspect influencing our re- 
sults is that, for the realistic band structure, the contributions 
from different segments of the Fermi surfaces to the net den- 
sity of states are weighted differently according to both the 
factor Nf :Tl in Eq. ( [23] ), and the segment length of the Fermi 
surface with a given direction of the Fermi velocity. Fig. [2] 



shows the profiles of the Fermi velocity and the correspond- 
ing weighting factors. The key point is that, in a tetragonal 
system, Vf )U can have a fourfold anisotropy in the plane that 
either enhances or competes with the gap anisotropy in deter- 
mining the contribution to the net DOS in the superconducting 
state, see Fig.[3ja). The detailed interplay of the two depends 
not only on the value of the angle-resolved DOS, but also on 
the direction of the Fermi velocity. 

Indeed, naively one might expect that the relatively large 
contribution to Nf^ n from the near-45° direction, combined 
with the node of the d x 2_ y 2 at the same angle in Fig. j3ja) 
should enhance the field-angle anisotropy for that symmetry 
of the superconducting state relative to the d xy case when the 
direction of the greatest Nf iU is fully gapped. In fact, the 
opposite is true, see Fig. [3jb): the oscillations are enhanced 
for d xy symmetry. 

This is an indication that the flat parts of the Fermi surface 
with large values of the Fermi velocity, see Fig. [2] contribute 
more to the total DOS, when the field is at 45° and all four 
flat parts are 'active'. When the field is along 0° or 90°, only 
two flat parts contribute. In contrast, the four 'active' corners 
with smaller velocities, and hence slightly larger Nf^ n (kf) 
give a smaller contribution simply because their arc length is 
a smaller fraction of the total Fermi surface length in the re- 
spective k z slice. It follows that d xy -pairing, which has nodes 
in the flat parts of the FS, exhibits enhanced C(a = 45°). In 
contrast, the d x i_ y i profile, gaps those regions of the Fermi 
surface and thus anisotropy of C is suppressed. Hence the ex- 
act role of the Fermi surface shape and curvature in the field- 
angle oscillations is highly non-trivial. 

At higher temperatures the simple low-T expression in 
Eq. (7) is only qualitatively correct, and both detailed 
calculation J 3 * 8 ! 42 ! 43 1 and experiment j 6 * 44 * 45 ' demonstrated that 
the anisotropy in the heat capacity is reversed relative to the 
low-T result. The lower panel in Fig. [3] shows the field- 
induced SC DOS as a function of quasiparticle energy below 
the SC gap for a = 0° and a = 45° for all three pairing 
symmetries considered here. We immediately see that the SC 
DOS at these angles switch and reverse magnitude, which re- 
flects in the sign reversal of the oscillations in specific heat as 
a function of temperature. Note that, due to the presence of 
Fermi velocities in k in Eq. (7), a one-to-one correspondence 
between SC DOS and k is not straightforward for FS that lack 
continuous rotational symmetry in the plane. 

We show below that for realistic and material- specific 
anisotropic FS, we still find the sign reversal of the heat ca- 
pacity oscillations for <i-wave pairing, which was previously 
reported for the rotationally symmetric cases. Hence this sign 
change is a generic feature of of nodal gaps. However, the key 
finding in this work is that for moderately anisotropic FSs, 
measurably large field-angle dependence in the heat capac- 
ity and thermal conductivity is obtained at high fields already 
for isotropic gaps, which can lead to misinterpretations if an- 
alyzed solely in this field range and in terms of simple har- 
monics of the SC pairing symmetries. We stress that multi- 
band effects add additional complexity to any analysis, due to 
competing FS anisotropies. For example, we have previously 
shown that if the FS anisotropy in different bands is opposite 
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FIG. 3. (color online) Fermi surface anisotropy of the normal- 
state DOS and SC gaps contrasted with the field-angle anisotropy 
of the Sommerfeld coefficient and the SC DOS. (a) The calculated 
FS anisotropy of the normal- state DOS juxtaposed with gap func- 
tions of three pairing symmetries. All the SC gaps are computed 
at the FS and all curves are shifted by their corresponding mini- 
mum value, except for the s-wave gap. (b) Specific heat coeffi- 
cient 7(a) = C(a)/T, normalized to its value at T c , calculated at 
T/Tco = 0.1 and H/H c2 = 0.1 for d-wave gaps and H/H c2 = 0.5 
for the s-wave gap. (c)-(e) Field-induced total SC DOS at T = vs. 
energy at two representative field angles a = 0° and 45° for all three 
pairing symmetries. Here we used H/H C 2 = 0.5 for s wave and 0.1 
for both d waves. Note the low- and high-energy crossings in the SC 
DOS (arrows) are related to the low- and high-T sign reversals in the 
oscillations of 7 and k in Fig. [4] 



(out-of-phase) to each other, then it can lead to additional sign 
reversals in the field-angle dependence of the thermodynamic 
quantities for s-wave gap, very similar to what was earlier ob- 
tained for nodal gaps only. 11 Furthermore, bands with differ- 
ent DOS lead to different amplitudes and shapes of the self- 
consistent value of the SC gaps (i.e., generally Ai 7^ A 2 ). In 
this case, the obtained numerical results become less intuitive 
to interpret and a simple one-to-one mapping between oscilla- 
tions and nodes is lost. However, at low temperature and low 
field the generic understanding of the anisotropy as a conse- 
quence of the nodal structure alone, remains valid. We give a 
detailed comparison of the different regimes below. 



B. Temperature evolution of field-angle-resolved oscillations 

We present the full angle-dependent profiles of 7(a) = 
C(a)/T and K,(a)/T for several temperatures at a representa- 
tive low field (H/H C 2 =0.1) for the two nodal gaps, and, for 
comparison, for an isotropic s-wave gap at a moderate field 



(H/H c2 = 0.5) in Fig. [5] In accordance with our earlier cal- 
culation for K y Fe2-xSz2 in Ref.lTTI we find that a substantial 
oscillation in 7 and k is present for isotropic s-wave pairing. 
The amplitude of the oscillation increases with stronger k z - 
dispersion. As in simple models, 8 close to the inversion line 
the oscillations are not a simple sum of the twofold and four- 
fold harmonics, but have a more complex profile. 

For nodal d xy and d x i_ y i pairings, the behavior of os- 
cillations of 7(a) and k{ol) is similar to results obtained 
for quasicylindrical FSs, 23 however the amplitude of oscil- 
lations and, crucially, the location of sign reversals in the 
H-T phase diagram are modified. Earlier such sign-reversal 
features were discussed o nly for highly anisotropic or nodal 
gap structures j 3 l 6 l 8 l 23 l 44 l 4 H Q ur material- specific results caution 
against straightforward interpretation of oscillations at inter- 
mediate fields as evidence of nodes, emphasizing the need to 
probe low energy excitations. 

We extract the amplitudes of the fourfold oscillations by 
defining 



where n£ = [C(a,T)/T]/[C N /T c ] and 

«4a(r) = p5 + n5o]/2- 



n 



45? 



(24) 



(25) 



where II£ = [n xx (a, T)/T]/[^/T c ], and C N and n N are 
the corresponding normal-state values at T c . This defini- 
tion removes any twofold, sixfold, etc., contribution from 
ft (a) originating from the field parallel or perpendicular to 
the vortex lines. 46 In fact, it is straightforward to show that 
for any function f(a) = Y^^=o a 2n cos(2na) the definition 



in Eq. ( |25) projects out any other harmonic contribution up 
to M = 5, resulting in = 2a^. We verified numerically 
that the amplitudes of twelvefold and higher order harmonics 



are negligible. On the other hand, the definition in Eq. ( [24] ) is 
less robust, but very convenient. It gives = 2a^ when 
a 2 = a 6 = aio = 0, which is sufficient when sample mis- 
alignment is negligible and when used away from the sign 
reversal line. 



C. H-T phase diagram 

In figure [5] we compile our results of the thermal quanti- 
ties into a contour map of the amplitude of the fourfold os- 
cillations in the H-T phase diagram for 7 = C/T (top row) 
and k/T (bottom row) for one nodeless and two nodal gaps. 
Recall that for quasicylindrical (rotationally-invariant in the 
basal plane) FSs the specific heat oscillations simply change 
sign between the d xy and d x i_ y i symmetries. 8 While the 
overall characteristics of the phase diagram remains qualita- 
tively the same for material-specific cases, substantial quanti- 
tative changes result from the inclusion of realistic Fermi sur- 
faces and the directional- and band-dependent contributions 
to the DOS. Important for the comparison with experiment, 
we find that the location of the sign-reversal lines for nodal 
gaps shown in Figs.|5ja2) and[5ja3) shifted compared to the 
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FIG. 4. (color online) Calculated oscillations of the heat capacity and thermal conductivity as a function of the field direction relative to the 
x axis, (al) Sommerfeld coefficient 7 = C/T normalized to its normal-state value Cn /T c q at fixed field H/H C 2 — 0.5 for s wave, plotted 
from low to high T in units of T/T c0 (bottom to top curves). Each curve is colored by the sign of the fourfold oscillation; a uniform color 
map is used for values below —0.025 and above 0.025. (bl) Same as in (al) but for the normalized thermal conductivity coefficient k/T. (cl) 
The fourfold amplitudes of 7 (dashed line) and k/T (solid line) are plotted as a function of T. For direct comparison the results for nodal d xy 
(H/H C 2 = 0.1) and d x i_ y i (H/H C 2 = 0.33) symmetries are plotted in panels (a2)-(a3), (b2)-(b3), and (c2)-(c3), respectively. Note that the 
non-vanishing of kau as temperature approaches the phase transition line in panels (cl) and (c3) is a consequence of the in-plane anisotropy 
ofH c2 . 



earlier simple models, due to the interplay of the SC order pa- 
rameter with the FS anisotropies. As a consequence, the sign 
of the fourfold oscillations, C^ a and ^4^, may be different 
over a wider range of temperatures and fields. This is to be 
contrasted with the results for rotationally symmetric Fermi 
surfaces, where the two were found to switch sign almost at 
the same temperatures and fields. We also verified that for s- 
wave pairing the high-T sign reversal is robust and remains 
at nearly the same location for a single-band superconductor 
with identical FS. 

Note also that at intermediate to high temperatures and 
fields there is very little in the heat capacity oscillation profile 



that distinguishes the isotropic gap from that of the d x 2_ y 2 
symmetry, see Fig. [5] panels (al) vs. (a3). However, there is 
a much more significant difference in transport, Fig. [5] pan- 
els (bl) vs. (b3), which implies that a simultaneous study of 
both C(a) and k(ol) is highly desirable to gain confidence 
about the underlying pairing symmetry in any multiband sys- 
tem where the low-temperature, low-field regime is experi- 
mentally unreachable. Of course, once the low energy sec- 
tor at low T and low H is accessed, the differences between 
different symmetries, and especially between the nodal and 
isotropic gaps, becomes obvious. Therefore, in general a 
rather detailed comparison between measurements and cal- 
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FIG. 5. (color online) Contour maps of fourfold amplitude oscillations of normalized specific heat C\ a (row a) and normalized thermal 
conductivity K^a (row b) in the H-T phase diagram. Each column denotes a different gap symmetry studied (isotropic nodeless s, nodal d xy 
and d x 2_ y 2 gaps). All plots use the same color map scale from minimum (red) to maximum (blue). Note that the fourfold amplitude is given 
with respect to field H || (100), i.e., a negative value corresponds to a minimum at a = 0°. Here T C {H) is defined by the vanishing of both 
gaps for given symmetry, which determines the line of the (minimum) upper critical field. Since the BPT approximation for isotropic s-wave 
pairing is not valid at low H, we shaded the corresponding area where our approach is not applicable. 



culations of the C^a and phase diagrams should be em- 
ployed to draw conclusions about the pairing symmetries. 



D. Comparison with experiments 

The superconducting Ce-1 15 compounds are well suited for 
the study of field- angle oscillations. Accordingly, there have 
been a number of experiments investigating the anisotropy of 
the thermal conductivity and the heat capacity under the ro- 
tated field. Here we compare the experimental results with 
our findings, previously shown in Fig. [5] Since the upper crit- 
ical field is Pauli-limited, and our calculation does not account 
for the Zeeman splitting, we cannot expect our results to map 
directly onto the measurements near H C 2. Nevertheless we 
believe that a qualitative comparison can be made, especially 
for systems with strong paramagnetism in the low-field part 
of the phase diagram, which, when rescaled to the appropriate 
values of the upper critical field, is essentially identical to that 
computed in the absence of the Zeeman term!^ 

Ce Coins : The unconventional nature of superconductiv- 
ity was recognized early on through the discovery of power- 



law dependence in the temperature behavior of the spe- 
cific heat and thermal conductivity!^^ magnetic penetr ation 
depth J2ED anc [ spin-lattice and muon-spin relaxation rate d 57 ! 58 ! 
consistent with predictions for nodal lines in the gap. On sym- 
metry grounds the anisotropy of the upper critical field van- 
ishes near T c0 , ff C 2(0°) = i^ C 2(45°). In our calculations, we 
find that noticeable anisotropy emerges for T/T c o < 0.7. In 
this range i^ C 2(0°) > H c2 (4:5 ) for both s and d x i_ y i pair- 
ing symmetries, while the anisotropy is opposite for d xy pair- 
ing, i.e., the nodal directions have a lower H c2 value. The 
anisotropy for s and d x i_ y i gap is in qualitative agreement 
with the H C 2 measurements of CeCoIn 5 by Settai et al.p^who 
reported i^ C 2(0°) > i^ C 2(45°) at low temperatures. The mea- 
sured anisotropy is only a few percent, which would be consis- 
tent with the assumption that the band electron ^-factor, and 
hence the Pauli limiting field is isotropic in the plane, and 
the weak anisotropy is due to a residual orbital effect. Re- 
markably, the opposite relationship, H C 2(0°) = 11.8 T < 
# C 2(45°) = 11.9 T was found in Ref.HI So far the experi- 
mental discrepancy of the in-plane H C 2 anisotropy remains an 
open puzzle. 

The original interpretations of the field-angle-resolved ther- 



10 



1) CeCoIn 5 ( a 2) CeRhIn 5 ( a 3) Celrln 5 




FIG. 6. (color online) Theoretical contour maps of fourfold amplitude oscillations for d x 2_ y 2 gap for specific heat (upper panel), and thermal 
conductivity (lower panel) of CeCoIns, reproduced from Figs.|5ja3)-b(3). The same theoretical data is repeated in three columns, but compared 
with three different experimental data for isostructural superconductors within the Ce-115 family, (al): Specific heat data of CeCoIn 5 by An 
et alP (circle s) an d Aoki et alP^ (squares), (bl): Thermal conductivity data of CeCoIn 5 by Izawa et alP^(a2): Specific heat data for CeRhIn 5 
by Park et alP^(a3): Specific heat data for Celrlns by Lu et al. 45 (circles) and Kittaka et alP (squares). (b3): Thermal conductivity data of 
Celrlns by Kashara et al. 49 The symbol size gives the corresponding reported amplitude of the oscillation, whereas the filled color depicts its 
sign. We find reasonable agreement between theory and experiment in both sign and amplitude of oscillations. 



mal conductivity 12 and specific heaP^ measurements contra- 
dicted each other regarding the location of the d-wave nodal 
lines in CeCoIn 5 . The controversy was finally settled by the 
observation of the inversion in the specific heat oscillations by 
An et alP In Fig.^al) and (bl) we plot both the C± a and 
experimental data points (symbols). The agreement between 
theory and experiment is quite convincing for 2 -wave 
symmetry and rules out pairing scenarios of either s or d xy 
gap. 

CeRhln^: The high-pressure, angle-resolved specific heat 
measurements of CeRhIn 5 by Park et alP^ showed a clearly 
delineated fourfold oscillation with C(0°) < C(45°), which 
was interpreted in favor of d-wave symmetry. The measure- 
ments were performed down to temperatures as low as 0.3 
K (0.3T/T C ) and in fields between 0.2 and 0.9 T. At the pres- 
sure of 1 .47 GPa the superconductivity coexists with antiferro- 
magnetism with superconducting transition T c = 1.04 K and 
in-plane H c2 = 1.2 T at 0.3 K. The measured in-plane H c2 
anisotropy was negligible. As we noted before, in this region 



of the H-T phase diagram both s-wave and d x i _ y i -wave gaps 
are nearly indistinguishable giving rise to fourfold oscillations 
with the minimum of C(a) occurring at H || (100). Support- 
ing the (i-wave interpretation, T- and i^-dependent measure- 
ments down to 0.3 K and 0.05 T exhibited power-law behavior 
consistent with unconventional superconductivity with nodes, 
i.e., C/T ~ T and C/T ~ y/H, respectively. Additional ex- 
periments at higher pressure (2.3 GPa), i.e., in the purely su- 
perconducting phase, and at T/T c = 0.14 and H/H c2 = 0.09 
showed evidence of fourfold oscillations with a negative am- 
plitude C^a of order 4%. 48 However, to unequivocally rule 
out the possibility of s-wave pairing, based on field-angle- 
resolved measurements alone, experiments would have to be 
performed at temperatures significantly below T c /3, where 
the exponential T-dependence of the fully gapped excitation 
spectrum becomes visible. Power laws were also seen in other 
pressure measurements of the specific heat, spin-l attice a nd 
muon-spin relaxation rates down to T/T c w 0.15! 58 * 61 ^ In 
chemically doped CeRhi-^Ir^Ins a T 3 dependence was seen 
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in 1/Ti just below T c , which tends toward linear in T at lower 
temperatures as is typical of dirty d-wave superconductors P^l 
In Fig.[6ja2) we plot the C^ a experimental data points (sym- 
bols) on top of the phase diagram for d x i_ y i gap. The field- 
angle-dependent experiments taken by themselves are incon- 
clusive, though combined with the reported T and H depen- 
dences are strongly suggestive of d x 2_ y 2 -wave superconduc- 
tivity in CeRhlns . 



Celrln^: There is an ongoing controversy about the pair- 
ing symmetry in this compound, because of its two differ- 
ent superconducting domes, namely one as a function of Rh 
doping and the other as a function of pressure. In addi- 
tion, there is disagreement over the interpretation of the ther- 
mal conductivity data. On one side, the field-angle-resolved 
measurements 49 and power-law dependence in temperature 
were argued as evidence for d-wave gap with vertical line 
nodes, similar to the sister compound CeCoIn 5 pl3 while on 
the other side thermal conductivity measurements, in partic- 
ular the temperature and magnetic field dependence of the 
residual value of k/T along different axes, were interpreted 
in favor of a three-dimensional hybrid gap with a horizontal 
line node The hybrid gap proposal was inspired by simi- 
lar gap functions studied some time ago for the heavy-fermion 
superconductor UPt3p2l22EH t f ur ther complicate the inter- 
pretation, the results by Shakeripour et al. were also argued 
to be consistent with vertical line nodesP^In addition, power 
laws were reported for magnetic penetration depth and spin- 
lattice-relaxation rate! 57 ! 63 ! 70 ' The temperature behavior of the 
anisotropic penetration depth was interpreted to be consistent 
with vertical line nodes but not with point nodes and a hor- 
izontal line node of the hybrid gap. 70 However, the conclu- 
sive evidence for the in-plane gap variation comes from very 
recent angle-resolved specific heat measurements at ambient 
and finite pressure. Lu et al. 45 (circles) reported fourfold os- 
cillations inside the pressure dome of Celrlns with sign rever- 
sal of the oscillations at high temperatures between 0.4 and 
0.6T C . These data taken together with a low-T anisotropy 
of H c2 (0°) > i7 c2 (45°) and the fact that this compound be- 
longs to the same family of Ce-115s was strongly suggestive 
of two-dimensional d x i _ y 2 -wave pairing with vertical line 
nodes. Unfortunately, the temperature in Ref.|45]was too high 
to formally exclude isotropic s-wave pairing, see the phase 
diagram in Fig. [5jal) versus (a3). The specific heat data of 
Ref. [71 on the other hand, were taken down to 80 mK, that 
is 0.2T C . Therefore, the specific heat oscillations are support- 
ive of the d x 2_ y 2 gap scenario. Data from both experiments 
are included in the comparison in Fig. [6ja3). In addition 
field-angle-resolved thermal conductivity data were reported 
by Kasahara et al.,E2 which are shown in Fig. (6^b3). Com- 
bined with the specific heat oscillations, they provide strong 
support for this pairing symmetry. Hence at present the over- 
whelming majority of experiments supports the ^2 -wave 
superconductivity with vertical line nodes in Celrln 5 . 



V. DISCUSSION AND CONCLUSIONS 

We performed realistic model calculations of the field- 
angle-resolved specific heat and thermal conductivity using 
a tight-binding parametrization of the electronic structure 
within a two-band model of superconductivity, which is rel- 
evant for the Ce-115 heavy fermions. Our systematic analy- 
sis of field-angle dependence showed that modest anisotropies 
in the density of states and the in-plane Fermi velocities of a 
tetragonal crystal contributes significantly to the fourfold os- 
cillations in the vortex state, when the magnetic field is rotated 
in the basal plane. As evidence we showed that such oscilla- 
tions exist at intermediate to high fields even for an isotropic 
5-wave gap. Remarkably, the sign reversal of fourfold oscil- 
lations occurs not only for nodal d-wave gaps, but also for an 
isotropic s-wave gap as the temperature is decreased. This is 
one of the main findings of this work and implies that away 
from the low temperature and low field region it may be dif- 
ficult to distinguish different pairing symmetries based on the 
field-anisotropy of a single probe alone. 

Finally, we compared the results of the field-angle-resolved 
calculations within our model with recent experimental data 
on different members of the Ce-115 family. The behav- 
ior of the self-consistently determined thermal quantities for 
nodal 6^2 -wave gap is consistent with experimental re- 
ports for CeCoIn 5 . The same phase diagram is also consis- 
tent with specific heat data for CeRhIn 5 and Celrln 5 . Since 
both CeRhlns and Celrlns have similar electronic structure 
as CeCoIn 5 near the Fermi energy, we believe that our Fermi 
surface parametrization is valid for all three compounds. Con- 
sequently, very similar phase diagrams should result for all 
three materials for which material- specific calculations dras- 
tically improved the agreement between theory and experi- 
ment. The comparison with experimental data is restricted 
to low fields , sinc e the superconductivity in this material is 
Pauli-limited,^ 4 -^ and there are indications of a quantum crit- 
ical point in the vicinity of the upper critical field at zero tem- 
perature, H C 2(0)W^ Thus the regime near the upper criti- 
cal field at low temperatures is beyond the scope of the cur- 
rent treatment. We find that within our realistic model of the 
Fermi surface parameters the fourfold anisotropy map is in 
better agreement with experiments on CeCoIn 5 , if we assume 
a weak dispersion along the k z axis. Note that the relatively 
small anisotropy of the upper critical field in this material does 
not have direct connection with the anisotropy of the electron 
dispersion, as it likely stems from the Pauli limiting of su- 
perconductivity. Since the electronic band structure is very 
similar among the Ce-115s near the Fermi energy, we expect 
that our findings for CeCoIns are also relevant for CeRhIn 5 
and Celrln 5 under pressure. Considering that questions re- 
main about the exact superconducting gap structure and po- 
tential spin-fluctuation nesting in the Ce-115s, 18 a definitive 
theoretical account of field-angle-resolved measurements is 
warranted. 

We conclude with a note of caution for interpreting field- 
angle-resolved oscillations. Our self-consistent two-band 
model calculations demonstrated that simple observations of 
oscillations and sign reversals in either C(a) or k{q) are not 
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direct evidence for the presence of nodes or minima in the 
gap structure. Such conclusions can be drawn either from 
low-energy measurements, or at higher temperatures and field 
from a comprehensive simultaneous analysis within the same 
framework of both C(a) and n(a) measurements. Only a sys- 
tematic analysis of the fourfold oscillations in the H-T phase 
diagram can constrain the space of possible pairing scenarios 
for a given material. 
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